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Abstract 

We investigated the influences of molecular crowding on biochemical reaction 
processes on two-dimensional surfaces, using the model of signal-transduction pro- 
cesses on biomembranes. We performed simulations of the two-dimensional cell- 
based model, which describes the reactions and diffusions of the receptors, signal- 
ing proteins, target proteins, and crowders, on the cell membrane. The signaling 
proteins are activated by receptors and induce target proteins to unbind from the 
membrane. We found that the reaction rates of two-dimensional systems consis- 
tently exhibit a maximum at a high volume fraction of molecules, such that two 
molecules in the vicinity cannot easily exchange their positions. We further demon- 
strated that molecular crowding influences the hierarchical molecular distributions 
throughout the reaction process. The signaling proteins tend to surround the re- 
ceptors, and the target proteins tend to become distributed around the signaling 
protein- receptor clusters. This distribution accelerates the receptor-signaling pro- 
tein and signaling protein-target protein reactions. Thus, molecular crowding fre- 
quently enhances reactions on two-dimensional surfaces, but restricts reactions in 
three-dimensional bulk systems. 
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INTRODUCTION 

Several living systems can sense and respond to environmental variations by means 
of internal biochemical processes. The adaptations of several cells and the cell fate 
of multicellular organisms (e.g., cell proliferation, differentiation, and apoptosis) 
are typical behaviors regulated by intracellular signal-transduction processes ([1]- 
[6]) . This simultaneous progression of several internal biochemical processes requires 
the synthesis and interaction of a number of different proteins on several biomem- 
branes and in the cytoplasm involving several macromolecules, the cytoskelton, and 
internal organs. 

Recent studies suggest that individual cells contain a much higher volume frac- 
tion of macromolecules than do cells in typical in vitro conditions ([71- I26j) . In fact, 
the total volume fraction of macromolecules in a typical cell is estimated to be 
50-400mg/mL. whereas that under typical in vitro conditions is estimated as 1— 
lOmg/mL (7). This situation of such a high volume fraction of molecules, com- 
monly called "molecular crowding", give rise to extreme spatial restrictions. Thus, 
the diffusion and deformations (reaction) of molecules in the cytoplasm are highly 
suppressed (f8 H2"Tj) . Such spatial restrictions are also expected to enhance protein 
folding (I22| I23p . formation and stabilization of intracellular architecture (24), and 
processive phosphorylation of erk map kinase (|25l |26|) . 

The transduction of signals from the extracellular environment starts with the 
activation of receptor and signaling proteins on the cell membrane. Thus, the 
sensing and response of cells is dependent on rapid and appropriate transport and 
reaction of signaling molecules in a two-dimensional space. Recently, imaging mea- 
surements of macromolecules on the cell membrane have been extensively performed 
(j27H37|) . In some of these measurements, subdiffusive motion was revealed to be the 
typical motion of membrane proteins (|34H37j) . This observation implies the exis- 
tence of some type of intrinsic membrane domain like raft or nonimaged molecules, 
which restrict the observed molecular motions by means of their excluded volumes. 
Thus, to uncover the functional performance and mechanisms of the upstream part 
of the signal-transduction processes, the influences of molecular crowding on the 
reaction and diffusion dynamics of two-dimensional systems should be investigated. 

Here, we investigated the influences of molecular crowding on biochemical pro- 
cesses, using an ideal model of typical signal-transduction processes on biomem- 
brane, which mimic the EGF-RAS-RAF and G-protein (guanine nucleotide-binding 
proteins) signaling process. This model consists of active and inactive receptors (R* 
and R), active and inactive signaling proteins (S* and S), binding and unbinding 
target proteins (T and T*), and nonreactive molecules (crowder, C), which diffuse 
and react in a two-dimensional space. The signal-transduction process is described 
as the cascade from the activation of receptors to the unbinding of activated target 
proteins from the membrane, through the following reactions (Fig. 1 ). 
(A) The active/inactive receptor autonomously becomes inactive/active with the 
reaction rate kn/kn*, respectively, described by 

R^R* (1) 



and 



R* ^R 



(2) 
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(B) When the inactive signaling protein makes contact with an active receptor, 
this signaling protein is activated with the reaction rate ks- The active signaling 
protein autonomously becomes inactive with the reaction rate kg. These processes 
are described by 

R + S % R + S* (3) 

and 

S* ^ S (4) 

(C) The target protein autonomously binds to the membrane with the probability 
Pi n . When a target protein makes contact with an active signaling protein, this 
target protein is activated with the reaction rate kx- Here, the target protein un- 
binds from the membrane as soon as it is activated. In the absence of activation, 
the target protein unbinds from the membrane with the reaction rate P ou t- These 
processes are described by 




S* +T — - ^> S* + T* (unbind from membrane) (6) 

and 

T owt ) T (unbind from membrane) (7) 



Here, each of molecules, R, R* , S* , S, T, and C has its own volume. These 
molecules move randomly under the restriction of their excluded volume, where 
the distance between the centers of two molecules cannot be smaller than the sum 
of their Van der Waals radii. For simplicity, we assume that the radius of each 
molecular is uniform. 

A recently recognized feature of molecular crowding is that the reaction activity 
of the system exhibits a maximum at an appropriate total volume fraction of the 
system (|9l H5 [T8]) . Such an appropriate volume fraction is estimated as not too 
large as to restrict molecular diffusion, but large enough to cause frequent collisions 
between reactive molecules. Recently, mean field analysis based on consideration 
of the balance between the diffusivity and collision rate of molecules has provided 
an explanation for the volume-fraction-dependent reaction activity of some three- 
dimensional bulk reaction systems ([9]). 

Here, we propose that the volume-fraction-dependent aspects of reaction pro- 
cesses on two-dimensional surface should deviate considerably from those of bulk 
reaction systems. Indeed, when analyzed by mean-field approximation, the present 
reaction system shows almost the same volume-fraction-dependent aspects as do 
those of the previous studies (fl~4"l [26]) (see Appendix). Here, the signal flow esti- 
mated as the rate of activation of T, exhibits a maximum value at an intermediate 
total volume fraction of the molecules. On the other hand, the simulation results 
of the present reaction system in a two-dimensional space indicate that the signal 
flow exhibits a maximum at a higher volume fraction of molecules, such that two 
molecules in the vicinity cannot easily exchange their position (see Fig. [4] Results 
and Discussion for details). 

The present simulation further show that, in the case exhibiting the maximum 
signal flow, spontaneous organization of spatially ordered molecular distribution 
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takes place. Thus, the signal proteins surrounding the receptors are rapidly ac- 
tivated, and the target proteins distributed around the receptor-signaling protein 
cluster and are rapidly activated. In this way, the formations of spatial structure 
that cannot be assumed in mean-field approximation provide the dominant con- 
tribution to the reaction activity. As mentioned in the following sections, such a 
spatial structure is formed and sustained by the typical processes of signal trans- 
ductions on the membrane, e.g., the binding, reactions, and unbinding of molecules, 
and also by the inferences of molecular crowding. 

MODEL AND SIMULATION METHOD 

Cell-based model 

Here, we have used the two-dimensional cell-based model ([551 ESI) to describe the 
diffusion and reactions of the active and inactive receptors (ill* and R), the active 
and inactive signaling proteins (S* and S), the binding target proteins (T), and 
the crowders (C) on the membrane. The space is divided into N two-dimensional 
hexagonal cells as shown in Fig. [2^[ Each cell can contain only one molecule, which 
represents the excluded volume effect among molecules. Each molecule randomly 
hops from one cell to a neighboring empty cell or reacts in the manner indicated by 
Eqs. [TJ [2j |4j [5] or [7J-|- Here, reactions [TJ [2j |4j and [7] occur spontaneously with each 
reaction rate, whereas two-body reactions [3] and [6] occur when two corresponding 
molecules exist in adjacent cells. In the empty cell, reaction [5] occurs on its rate 
Pi n ; here, the crowder does not react. 

As the index of molecular crowding, we define "occupancy" as [Total number of 
molecules in the system] /N. The volume fraction, which has frequently been used 
as the index of molecular crowding in recent arguments, is roughly estimated from 
the occupancy, under the following assumptions. First we assume that all molecules 
are almost the same size. Second, the size of each cell is considered sufficiently large 
compared with those of the molecules, that indicate, the center of each molecule 
can fluctuate in each cell, and two molecules in neighboring cells do not maintain 
contact with each other. Third, the size of each cell must be small enough for the 
distance between two molecules in neighboring cells to be always smaller than the 
molecular diameter, to avoid the possibility of invasion of a molecule between two 
neighboring molecules. For example, we assume that the length of the edge of each 
hexagonal cell is ~ 1.7x [molecular diameter], which is close to the largest value 
to satisfy the above conditions. By this assumption, the volume fraction of the 
most dense condition with occupancy = 1 in the two-dimensional cell-based model 
is estimated as ~ 41% which is slightly large, but expected as the possible value in 
experimental situations (jTj). 

If we regard the present model as the G-protein signal-transduction system, its 
time and spatial scales are estimated by considering the size of each cell to be of 
the same order as its molecular size, i.e., nm. The diffusion rate of particles was 
estimated as ~ 10 _1 -10 _2 um 2 /s = 10 _1 -10 _2 nm 2 /us for the G-protein-coupled 
receptor without the cytoskeleton (1340 . Hence, the time interval between which 
a molecule moves a single molecule length is estimated as ~ 10-100 \is. We also 
assume that the 0.01-1 reactions occur for each molecule at each time step. This 
assumption gives the characteristic time for reactions as 0.1-10ms which is consis- 



Enhancement of Signaling Processes 



5 



tent with the time scale of the conformational change of typical proteins. 

Simulation method 

To simulate the present cell-based model, we used the Monte-Carlo method. The 
temporal evolution of the system is progressed by the iteration of the following 
steps. 

(i) One of the cells is chosen randomly. 

(ii) If this cell contains a molecule R, R*, S*, or T, the corresponding reactions, [TJ 
[2j |4j or [7j occur with the respective reaction rates. Here, reaction [7] indicates that 
T is removed from this cell. 

(hi) If this cell contains a molecule S or T, the corresponding two-body reactions, 
[3] or |6j occur with the rate given by the product of its [reaction rates] and [the 
number density of the corresponding catalyst on (six) neighboring cells]. Here, T* 
is removed from this cell as soon as it appears. 

(iv) If this cell contains a molecule, but no reactions occur, this molecule moves to 
one of the neighboring (six) cells chosen randomly when the chosen cell is empty. 

(v) If this cell contains no molecules, the binding process of a target protein [5] 
occurs with the reaction rate Pi n . This indicates that T is bound to this cell with 
the reaction rate P{ n . 

Here, the unit time step is given by the [Number of above iterations] / N . We 
dehne the signal flow of the system, J, as the frequency of activation to unbind 
the target protein, where J at time t is given by [Number of activations of target 
proteins between t and t + 1] . 



RESULT AND DISCUSSION 

Simulation result 

In this section, the typical properties obtained through the simulation of the model, 
which the system includes no crowders, are considered. Firstly, we focus on the 
steady-state signal flow J, the frequency of activation of the target proteins, for 
several volume fractions of the signaling protein [Stot] = [S] + [S*] and the binding 
probability of the target protein Pj n . Here, [X] means the volume fraction of the 
molecule X in the system. For simplicity, some parameters are fixed as follows: 
k s = 0.3, k s * = 0.3, k T = 0.3, and [R tot ] = [R] + [R*] = 0.25. However, the 
following arguments are qualitatively independent of these details. 

Fig. [3 | shows J as a function of [Stot] and Pi n for the set (kn,kji*, P ou t) = 
(A) (1,0,0), (B) (0.5,0.5,0), (C) (0,0,0.3), and (D) (0.5,0.5,0.3), Fig. g]A gives 



the cross section of Fig 3 A which showing J as a function of Pi n for several [Stot] 
(0 < [Stot] < 0.75, because [S tot ] + [Rtot] < 1 and [R to t] = 0.25). Together, these 
figures indicate that the signal flow J exhibits a maximum at Pi n = 1 with the 
respective appropriate values of [Stot]- The cases of Pj„ = 1 are such that a target 
protein binds to the membrane as soon as an empty space appears. 

This further means that J exhibits a maximum in the case of a much higher 
total volume fraction of molecules, p, than that expected by mean-field analysis (see 
Appendix and Fig. [TT"|B). Here, at [Stot] = 0.25,0.35, and 0.65, J increases linearly 
with p and exhibits a maximum at the maximum p for a given parameter set. It 
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should be noted that the diffusion of each molecule on the membrane is highly 
suppressed because most of the space is occupied by molecules in such cases of p 
such that J exhibits a maximum. Indeed, in the case of small [Stat], i- e -> [Stot] = 
0.025, 0.1, J has a maximum peak at low p, such that each molecule can easily 
diffuse. However, influences of the molecular crowding, which generally suppress 
the reaction activity in three-dimensional bulk systems can lead to predominantly 
larger signal flows in the present reaction processes. 



Spatial organization 

As shown in above, the numerical results of the present model deviate considerably 
from those expected by mean- field analysis. This finding implies the existence of a 
spatially nonuniform distribution of molecular species, which can cause the higher 
signal flow enabled by a high volume fraction of molecules. Thus, to observe the 
characteristic spatial distributions of molecules, we measure the radial distribution 
function of each molecular species around each receptor. 

The radial distribution functions of the receptor, signaling protein, and target 
protein around the receptor, dn(r), ds(r), and dr{r), are defined as 

, / \ (pR(r)+PR* (r)) , , 

d n( r ) = F^l . (8) 

ds(r) = te(r ' + 7 (r)) , (9) 
Mr) = (10) 

respectively. Here, px{f) (X = {R, R* , S, S* , T}) indicates the volume fraction of 
molecules X in the region where the distance from the receptor is r, which is given 
as the minimum number of steps to move from the receptor to the molecule X. ( ) 
means the sample and the long-time averaged value. The molecule X is considered 
to accumulate when dx (r) > 1, whereas X is considered sparse when dx(r) < 1. 

Firstly, we focus on the radial distribution functions when the system realizes 
the maximum signal flow. Fig. [5[] shows the radial distributions of the receptor 
(plus), signaling protein (cross), and target protein (circle) as a function of distance 
from a receptor, r, with a large volume fraction of molecules such that J exhibits a 
maximum at (kji, , P ou t) = (1,0,0). This figure contains two regions divided by 
r = r c ~ 5. For r < r c , the distribution of target protein tends to be sparse, whereas 
receptors and signaling proteins tend to accumulate. On the other hand, for r > r c , 
all molecules are uniformly distributed. Thus the molecules tend to distribute 
according to the following spatial structure: the signaling proteins surround the 
receptors, and the target proteins become distributed around the receptor-signaling 



protein cluster (Fig. 6 ) 



If molecular distribution occurs according to the abovementioned structure, S 
around R* and also T around S*, tend to be activated rapidly. Thus, the reaction 
process of the present model progresses actively, even at a much higher volume 
fraction of molecules. 

Moreover, the progress of such reaction processes is considered to contribute to 
the formation and sustainment of radial distributions in the following way. Once a 
signaling protein neighbors a receptor, the two tend to remain in the vicinity with 
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each other because the volume fraction of molecules is too high to diffuse them. 
Thus, as soon as the receptor is activated to be transformed into R, the signaling 
protein is activated to be transformed into S. Next, a T around this S is rapidly 
activated and unbound from the membrane. Soon, another T comes close to this 
S and is activated because the volume fraction of molecules is too high. These 
processes induce the flow of molecules from the distant location to the vicinity of 
the activated R, thereby sustaining the spatially ordered molecular distribution. 

It should be noted that the formation of the present spatial structure is induced 
by the typical processes of signal transductions on biomembranes, e.g., binding, re- 
actions, and unbinding of molecules, and also by the influences of molecular crowd- 
ing. If the volume fraction of molecules is too small to allow smooth diffusion on 
the membrane, the abovementioned spatially ordered molecular distribution does 
not emerge. 

Fig. [7J shows the radial distributions of each molecular species at r = 1 as a 
function of Pi n . With decreases in Pj n , regarded as decreases in the total volume 
fraction of molecules, the radial distributions at r = 1 become deviated from the 
case of a larger Pj n . Although the target proteins tend to exist in the vicinity 
of the receptors, such a configuration does not contribute positively to signal flow 
because the receptors cannot directly activate the target protein. Moreover, such 
a distribution may often disturb the activation of signaling proteins by activated 
receptors. Thus, the reaction activity decreases in line with a decrease in Pj„ 
because of the breaking of the spatially ordered molecular distributions. 

In other words, molecular flow is enhanced not diffusion, but by the spatial 
molecular ordering caused by influences of the molecular crowding. Thus, the re- 
action activity is increased. 



Reaction system with crowing molecules 

In the previous sections, we considered the ideal system, which is assumed to com- 
prise only the components of the considered signal-transduction process. How- 
ever, in general, several reaction processes are simultaneously progressed by several 
macromolecules on the cell membrane. The components progressing such a reac- 
tion process often behave as obstacles for the components progressing alternative 
reaction processes. Thus, to elucidate the influences of molecular crowding using 
a more realistic model of signaling processes on the membrane, we performed the 
simulation of the system containing the crowder molecule C, which moves randomly 
on the membrane without the reaction, binding, and unbinding. Here, C hinders 
the random movements of other reactive molecules because of its excluded volume. 
We consider the case of [Rtot] = 0.25 and [C] = 0.3, which is larger than or similar 
to the order of [Stot]- 



Fig. 8 shows J as a function of [Stot] and Pin for the set (/cr, kp*, P ou t) = (A) 



(1,0,0), (B) (0.5,0.5,0), (C)(l, 0,0.3), and (D) (0.5,0.5,0.3). Fig. [9J A and B 
and Fig. |10 [ A and B indicate the cross sections of Fig. [8[j A and C for some 
[Stot] (0 < [Stot] < 0.4, because [S tot ] + [Rtot] + [C] < 1). For the case of small P out , 
as shown in Fig. 8 A and B and Fig. [9~|A, J has a maximum peak at the certain 
[Stot] and Pi n , which is much smaller than 1, in contrast to the system without C. 
However, Pj n and the total volume fraction of molecules that exhibit a maximum 
J value are still much higher than those expected by mean-field approximations 
(Fig. 1 1 1 | B). Moreover, in contrast to the results of mean- field approximations, J 
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continuously changes, and has a finite value is maintained at Pi n = 1. 

These data are explained as follows. As mentioned in the previous sections, 
the flow of molecules toward the receptor-signaling protein clusters (R-S clusters) 
emerges if the present reactions proceed in the presence of a high volume fraction 
of molecules. In such a case, not only the target proteins but also the crowders 
flow toward the R-S clusters. It should be noted that the target proteins tend to 
be unbound as soon as they come into contact with the R-S clusters, whereas the 
crowders tend to remain bound. With increase in the volume fraction, the crowders 
come close to the R-S clusters, thereby suppressing the activation of target proteins. 

On the other hand, if the target proteins are frequently unbound from the 
membrane without activation, the emergence of empty spaces often enhances the 
exchange of positions between the target proteins and crowders. Thus, the sur- 
rounding probability of R-S clusters by crowders decreases, and target proteins can 
come close to the R-S clusters to be activated. The unbinding of target proteins 
without activation inevitably results in a certain loss of signal flow. However, based 
on an appropriate unbinding rate of target proteins P ou t, we obtain the signal flow 
J, which exhibits the same qualitative properties as the J value obtained in the 



system without C (Fig. 3 C and D). It should be noted that it is natural for the 
model of reaction processes on biomembranes to assume that the unbinding rate 
of the target protein without activation has a finite value. Thus, molecular crowd- 
ing generally tends to enhance the reaction activity in two-dimensional reaction 
systems such as signaling processes on biomembranes. 



CONCLUSION 

We investigated influences of the molecular crowding on the activity of reaction pro- 
cesses on two-dimensional surfaces, using the cell-based model of signal-transduction 
processes on the biomembranes. By simulation of the model based on the diffusion 
and reaction among receptors, signaling proteins, target proteins, and crowders on 
two-dimensional surface, we found that the signal flow exhibits a maximum value 
at a high volume fraction of molecules such that the diffusion of molecules is highly 
suppressed. We further demonstrated that such molecular crowding influences the 
hierarchical molecular distributions throughout the reaction process. The signaling 
proteins tend to surround the receptors, and the target proteins and crowders tend 
to become distributed around the receptor-signaling protein clusters. This distri- 
bution accelerates the receptor-signaling protein and the signaling protein-target 
protein reactions. The formation of such a spatial structure is induced by the re- 
striction of diffusion by molecular crowding and also by the typical processes of 
signal transductions on biomembranes, e.g., binding, reactions, and unbinding of 
molecules. In this way, molecular crowding frequently enhances the reactions on 
two-dimensional surfaces, such as signaling processes on biomembranes, but sup- 
presses reactions in three-dimensional bulk system. 

To control the signal transduction activity on the membrane, we expect the 
formation of the presented hierarchical molecular distributions provide the dom- 
inant contribution as well as the receptor clustering (|4U|) . On the other hand, a 
large number of reaction processes other than the signaling cascade are known to 
progress on the interior and exterior surfaces of several biomembranes of internal 
organs such as the mitochondria, Golgi body, and nucleus. Based on our present 
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argument, we predict that the characteristics of such intracellular reaction dynam- 
ics in two-dimensional systems are effectively lower than those in three-dimensional 
space. 

ACKNOWLEDGEMENT 

The authors are grateful to Y. Togashi, K. Takahashi, K. Aoki, M. Nishikawa, M. 
Kikuchi for useful discussion and information. This work was supported by the 
Research Fellowship of the Japan Society for the Promotion of Science for Young 
Scientists to MF, by the Global COE Program G14 (Formation and Development 
of Mathematical Sciences Based on Modeling and Analysis) of MEXT of Japan to 
HN, and by the Grant-in- Aid for Scientific Research on Innovative Areas (Spying 
minority in biological phenomena (No. 3306) (24115515)) of MEXT of Japan to AA. 



APPENDIX: MEAN-FIELD APPROXIMATION 

We analyzed the present model by the mean-field approximation, in a similar way 
to recent studies (|14| [26]) . For simplicity, we showed only the results of the system 
without crowders because the results are qualitatively the same between systems 
with and without crowders. To compare the simulation results in the text, we 
used occupancy, and not volume fraction, as the index of molecular crowding. The 
relationship between occupancy and volume fraction is described in the text. The 
temporal evolution of the density of each molecular species, [R], [R*], [S*], [S] and 
[T], is obtained by the modified mass action law involving the influences of the 

-k R [R] + k R *[R*] (11) 

k R [R] - k R *[R*] (12) 

-k s f(pW][S\+ks.[S*] (13) 

ksf(p)[R*}[S}-k s *[S*] (14) 

f(p)P m - k T f( P )[T][S*] - P out [T] (15) 
[Stot] + [Riot] + [T}. (16) 

Here, f(p) denotes the occupancy, p, dependency of the collision rate of two 
reactive molecules. According to recent studies, f(p) is chosen as a monotone 
decreasing function of p and implements /(0) = 1 and /(l) = 0, to represent the 
crowding effect (fT4"l |2U|) . Here, we employ f(p) = 1 — p and consider the case of 
which the receptors are always active and the target proteins do not autonomously 
unbind from the membrane, i.e., [Rtot] = [R*] an d P ou t = 0. We can obtain the the 
same qualitative results independently of the details f(p) and parameter sets. 

From these equations, we estimate the steady-state signal flow, defined as the 
frequency of activation for unbinding of the target protein, J = kTf(p)[T][S*]. Fig. 
|11 [ A and B shows J as a function of the binding probability of the target protein, 



occupancy of molecules. 

d[R] _ 
~dT ~ 
d[R*} _ 
dt ~ 

4S]_ 

dt 
d[S*} _ 

dt ~ 
d[T] _ 

dt 

P = 
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Pin, and of the total volume fraction of molecules on the membrane p, respectively. 
In Fig. |11 | A, the maximum peaks in J for several P{ n are obtained in a very small 
region (~ 10 3 ), and each J drops discontinuously to zero at each critical value. 
Moreover, as shown in Fig. |11 | £>, there exists a suitable intermediate value of 
total volume fraction of the molecules and a suitable mixing rate of each molecular 
species, to maximize J. 
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FIGURE 1 : Signaling pathway of this study. The signal is transferred via the activation 
of the receptor (Eq. [T]), the activation of the signaling protein (Eq. [3]), and the activa- 
tion of the target protein with unbound (Eq. [6]). The target protein in the cytoplasm 
stochastically becomes inactive and is rebound to the membrane (Eq. |5k; the receptor 
and signaling protein also stochastically becomes inactive (Eqs. [2] and [4]). The target 
protein on the membrane is stochastically unbound from the membrane (Eq. [7]). 




FIGURE 2 : Outline of the study model. The membrane is described as a hexagonal 
lattice surface. One cell can contain only one protein. R, R* , S, S* , T, and C indicate 
the inactive receptor, active receptor, inactive signaling protein, active signaling protein, 
target protein, and crowder, respectively. 
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FIGURE 3 : Amplitude of signal flow J (color gradation) for (kn, k^*, P ou t) = (A) 
(1.0,0,0), (B) (0.5,0.5,0), (C)(1.0, 0,0.3), and (D)(0.5, 0.5, 0.3); and [R tot ] = 0.25 as a 
function of the volume fraction of the signaling protein [Stat] (horizontal axis) and the 
binding probability of the target protein P in (vertical axis), respectively. The brighter 
color indicates the higher signal flow. 
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FIGURE 4 : Cross section of Fig. 3 A., showing J as a function of (A) Pi n and (B) p for 
[Sua] = 0.025 (red, +), 0.1 (green, x), 0.25 (blue, □), 0.35 (purple, O), and 0.65 (yellow, 
A). 
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FIGURE 5 : Radial distributions around the receptor as functions of distance from 
receptor, r, for [R to t\ = 0.25, [S tot ] = 0.35, and P in = 1, the combination realizing a max- 
imum signal. The red (+), green (x), and blue (O) lines indicate the radial distribution 
of the receptor (d^(r)), signaling protein (ds(r)), and target protein (dr(r)) around the 
receptor, respectively. 




FIGURE 6 : Schematic diagram of spatial structure in the crowded condition with a 
high rate of signal flow. 
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FIGURE 7 : Radial distributions around the receptor for r = 1, d R (l), d s (l), and 0^(1), 
as functions of P in when [R] = 0.25 and [S tot \ = 0.35. The red, green, and blue lines 
indicate c(r(1), ds{l), and rfr(l), respectively. 
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FIGURE 8 : Amplitude of signal flow J (color gradation) for (&r, &_r*, P ou t) = {A) 
(1.0,0,0), (B) (0.5,0.5,0), (C)(1.0, 0,0.3), and (D)(0.5, 0.5,0.3); [R to t] = 0.25, and [C] = 
0.3 as a function of the volume fraction of the signaling protein [Stot] (horizontal axis) 
and the binding probability of the target protein P in (vertical axis). The brighter color 
indicates the higher signal flow. 




FIGURE 9 : Cross section of Fig. 8 A, showing J as a function of (A) Pi n and (B) p 
for [Sua] = 0-025 (red, +), 0.1 (green, x), 0.25 (blue, □), and 0.35 (purple, Q). 
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FIGURE 10 : Cross section of Fig. |8p, showing J as a function of (A) Pj n and (B) p 
for [Sua] = 0.025 (red, +), 0.1 (green, x), 0.25 (blue, □), and 0.35 (purple, CO- 




FIGURE 11 : Signal flow J obtained by mean-field approximation as a function of (A) 
P in (xl0~ 3 ) and (B) p. [S tot ] = 0.1 (red, +), 0.2 (green, x), 0.3 (blue, □), 0.4 (purple, 
O), and 0.5 (yellow, A). 



